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Abstract 

This article introduces the splitting method to systems responding to rough paths as ex- 
ternal stimuli. The focus is on nonlinear partial differential equations with rough noise but we 
also cover rough differential equations. Applications to stochastic partial differential equations 
arising in control theory and nonlinear filtering are given. 

1 Introduction 



This article introduces the splitting-up method for systems responding to rough paths as external 
stimuli. We deal with (nonlinear) rough partial differential equations, RPDEs, formally written as 

du = F (t, x, u, Du, D 2 u) dt + A (t, x, u, Du) dz on (0, T] x R e , u (0, x) = u Q (x) (1) 

but we also cover rough differential equations, RDEs, of the form 

dy t = V(y t )dt + W{y t )dz t . 

In both examples, z is an external stimuli given as a rough path, F is a nonlinear (possibly degen- 
erate) elliptic operator, A is a collection of affine linear operators, i.e. 

A fe (t,x,r,p) = (p-a k {t,x)) + rv k (t,x) + g k (t,x) . 

and a, v, g resp. V, W = (Wj) are (collections) of vector fields on [0, T] x R e resp. R e . Consequences 
of our results are splitting results for (nonlinear) stochastic partial differential equations, SPDEs, 
that is, when z is taken to be a rough path lift of a stochastic process (e.g. fractional BM with 
Hurstparameter > i); for example, linear SPDEs with Brownian noise in Stratonovich form (e.g. 
the Zakai equation from nonlinear filtering), 

d 

du = L (t, x, u, Du, D 2 uj dt + A k (t, x, u, Du)odB k , u(0,.) = M (■) • (2) 

fc=i 
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are covered. L is here a linear (degenerate elliptic) operator, 



L (t, x, r, p, X) = Trace [A (t, x) ■ X] + b (t, x) ■ p + c (t, x, r) , 

and B a standard G?-dimensional Brownian motion. The splitting-up method (which runs under 
many names: dimensional splitting, operator splitting, Lie-Trotter-Kato formula, Baker-Campbell- 
Hausdorff formula, Chernoff formula, leapfrog method, predictor-corrector method, etc.) is one of 
the most prominent methods for calculating solutions of (stochastic-, ordinary-, partial-) differential 
equations numerically; for a survey we recommend |MQ02 . For S(P)DEs a splitting-up method 
was introduced in [BG89] for the Zakai equation in filtering and has received much attention since. 
We explicitly mention |KG03j which extends the previous results to general linear SPDEs of the 
form ©. 

Let us informally describe the general idea of splitting using ([2]) : For neN consider the partition 
D n = {<" = in~ 1 T, i = 0, . . . , n\ of the interval [0, T] and define the approximation u n recursively 

U n (C+l, ■) := [Qt« t?+1 o P t?t?+1 ] (U n (t?, ■)) 

with {P s t, < s < t < T} and {Q s t, < s < t < T} the solution operators of 

dv — L (t, x,v, Dv, D 2 v) dt, v(s,x) — v(x) and (3) 

d 

dw = Afc (t, x, w, Dw, D 2 w) o dBj, w (s, x)=v (x) . (4) 

k=l 

That is, on each interval one solves first equation ([JJ on with initial value 

u n (t™, .) and then one uses its solution as initial value for the PDE ([3]) (so-called "predictor" and 
"corrector" steps in [FLG91] ). Under appropiate conditions, one can show that u n converges to u 
and also derive rates of convergence, |KG03j . 

All the above mentioned authors use (to the best of our knowledge) either semigroup theory or 
stochastic calculus to prove splitting results but neither are available for ([T]) due to the nonlinear 
operator F and the nonsemimartingale noise z. The point of view of this article is different; loosely 
speaking: splitting-up results follow from stability in a rough path sense. We combine the method 
of Krylov and Gyongy, [KG03], of stretching out the time-scale with certain stability results of 
RPDEs. Applications to SPDEs then follow and our results are, to the best of our knowledge, new 
for nonlinear PDEs with noise of above form (see also |LS98al ILS98bl lLS98bl ILSOOal ILSOObj ) . Due 
to the generality of equation {l} we do not give rates of convergence but hope to return to this 
question in the future. 



1.1 Some ideas from rough path and viscosity theory 

Let us recall some basic ideas of (second order) viscosity theory [CIL921 IFS06bj and rough path 
theory LQ02] ILGL07j . As for viscosity theory, consider a real- valued function u — u(t,x) with 
t e [0, T] , x £ M. e and assume u 6 C 2 is a classical subsolution, 

d t u + F (t, x, u, Du, D 2 u) < 0, 

where F is a (continuous) function, degenerate elliptic in the sense that F (t, x, r, p, A + B) < 
F (t,x,r,p, A) whenever B > in the sense of symmetric matrices. The idea is to consider a 
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(smooth) test function ip and look at a local maxima (t, xj of u — ip. Basic calculus implies that 
Du (t, xj = Dip (t, xj , D 2 u (t, x^j < Dip (t, x^j and, from degenerate ellipticity, 

d t ip + F (t, x, u, Dip, D 2 cp) < 0. (5) 

This suggests to define a viscosity supersolution (at the point (x, tj) to d t + F = as a continuous 
function u with the property that (JS|) holds for any test function. Similarly, viscosity subsolutions 
are defined by reversing inequality in ([5]); viscosity solutions are both super- and subsolutions. A 
different point of view is to note that u (t, x) < u (i, xj — p (i, xj + ip (t, x) for (t, x) near (t,i). A 
simple Taylor expansion then implies 

u (t, xj < u (t, xj + a (t — tj + p ■ (x — xj + — (x — xj T ■ X ■ (x — x) + o (\x — x\ 2 + \i — t^j (6) 

as \x — x\ 2 + \i — 1 1 — > with a = dtp (i,x) , p — Dip (t, x) , X = D 2 p (i, x) . Moreover, if ([5]) holds 
for some (a,p, X) and u is differentiable, then a = dtu (i, x) , p = Du (t, x) , X < D 2 u (t, x), hence 
by degenerate ellipticity 

d t tp + F(i,x,u,p,X) < 0. 

Pushing this idea further leads to a definition of viscosity solutions based on a generalized notion 
of "(dtu, Du, D 2 u) v for nondifferentiable u, the so-called parabolic semijets, and it is a simple 
exercise to show that both definitions are equivalent. The resulting theory (existence, uniqueness, 
stability, ...) is without doubt one of the most important recent developments in the field of partial 
differential equations. As a typical result, the initial value problem {dt + F)u — 0, u (0, •) = 
u e BUC(R e ) has a unique solution in BUC([0,T] xl e ) provided F = F{t,x,u,Du,D 2 u) is 
continuous, degenerate elliptic, proper (i.e. increasing in the u variable) and satisfies a (well-known) 
technical conditional. In fact, uniqueness follows from a stronger property known as comparison: 
assume u (resp. v) is a supersolution (resp. subsolution) and uq > vq; then u>»on [0,T] x R e . 
A key feature of viscosity theory is what workers in the field simply call stability properties. For 
instance, it is relatively straight-forward to study (dt + F) u = via a sequence of approximate 
problems, say (dt + F n ) u n — 0, provided F n — > F locally uniformly and some apriori information 
on the u n (e.g. locally uniform convergence, or locally uniform boundednesqj- Note the stark 
contrast to the classical theory where one has to control the actual derivatives of u n . 

The notion of stability is also central to rough path theory. Let yo S M e , V, W = (Wi) i=1 d be 
(collections of) vector fields on R e . Using rough path theorjQ, one can speak of solutions to 

dy t = V (y t ) d£ + W (y t ) dz t , y (0) = y G M e (7) 

if the weak geometric rough path z gC""™ ([0, T] , (R d )) and f G C 9 " var ([0, T] , K) have "com- 
plementary Young regularity" - + - > 1: then (£, z) can be seen as "time-space" rough path, 

i.e. an element of C max (p>«) ([Q, T] , G[ max (p>«)] (r © R d ) \ since all necessary cross iterated integrals 
between z and £ are well-defined (using Young integration; [FV09 , LV06]). Now any sequence 

1 BUC (...) denotes the space of bounded, uniformly continuous functions. 

2 (3. 14) of the User's Guide [CIL92] , 

3 What we have in mind here is the Barles-Perthame method of semi-relaxed limits FS06b . 
4 Young theory is actually good enough for p < 2; see |Lyo94| . 
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(C,z n ) n C C 1 "™ 1 ' ([0,T] ,R) x C 1 "™ 1 ([0,T] ,R d ) leads to a sequence of ODE solutions (y n ) n of 

(again assuming the vector fields are regular enough) and we call any limit point in uniform topology 
of (y n ) n a solution of the RDE 

dvt = V(vt)dZ t + W(yt)dxt, 
if (£", z n ) converges to (£, z) in the sense 

^P||5 [p] (z")|| p _ var + sup||5 [9] (r)||„ < oo (8) 

n n 

do (S [p] (z n ) , z) + \C - CL;[o,T] -> as n ^ oo. 

Here, 5r p i denotes the canonical lift to a geometric step-[p] rough path (i.e. given by Riemann- 
Stieltjes integration) , 

S]p] {z n ) s ,t = 1 + f d< +•••+/ rf< ® • • • ® < w ■ (9) 

«/ S Js<1ti<--<l[[ p ]<f 

If there exists a unique solution we denote it by nv,(W) {^i x \ (^j z )) to emphasize dependence on 
the initial condition yo, the rough path (£, z), the vector field V and the collection of vector fields 
W = (Wi)f =1 . A variation of Lyons' limit theorem, |FV09j . gives sufficienl(3 conditions for existence 
and uniqueness of such RDE solutions. 

2 Two examples 

In this section we sketch our approach on two examples, ODEs and RDEs (resp. SDEs when take 
the rough path lift of a stochastic process). Therefore let us introduce some notation: for fixed 
A > and t £ [0, T] sei0 t A = [t/A\ A and t A = [t/A\ A + A (i.e. [t A ,t A ) is th e interval in the 
partition of [0, T] of constant mesh size A that contains t). Motivated by |KG03j define two time 
changes 

/ a j.\ / t A +2(t-t A ), t A <t<t A + A/2 / A\ 

a(A ' t) = { t A + A/2<t<t* ,H^t)=a(A,t+-y 

That is, a (A, .) runs on the first half of each intervall [t A ,t A ) with double speed from t A to 
t A and stays still in the second half, wheras b (A, .) does this in opposite order; also the paths 
(a (A, ■) , b (A, •)) converge in (1 + e)-variation against the path icfe : t i — > (t, t) for every e > 0. 

Further, for given n > 1 denote by D n the partition {%T, k = 0, . . . , n) of [0, T]. 

The following two examples were shown to us by Terry Lyons at the IRTG SMCP Summer 
School 2009 in Chorin: 

5 ...and essentially sharp (see IDav07n . We remark also that one can slightly improve on regularity assumption 
of the vector fields by reproving the Universal Limit theorem for RDEs with drift instead of using the "space-time" 
rough path (z, £) 

e [.J denotes the lower floor function. 
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2.1 Splitting ODEs 

Let V, W g Lip 1 (R e ,R e ). We are interested in splitting of the ODE 

dy t = V (y t ) dt + W (y t ) dt, y (0) = y G R e . (10) 

As in the introduction, denote the solution of (|10[) as ttv.w (OjJ/o!*^)- Classic Lie splitting corre- 
sponds to the approximation of the path idi by the sequence of paths 1 1 — > (a t) , 6 i)) . 
Therefore let y™ be the ODE solution of 

dyt =V(y t )da(n- 1 ,t)+W(y t )db(n- 1 ,t), y (0) = y g R e , 

i.e.?/" = 7rv,w (0, yol (a(re _1 , •) , 6 (n _1 , •))). For brevity, define the solution operator |p™' t v ,0 < s < £ < t| 
and {Pjf , < m < T} mapping points in R e to R e as 

P™'^ (a;) := 7ry (s, x; a ■)) and Pjl s (ac) := 7ry (0, x; id\) t _ s = 7ry (s, x; id\) t 

(here idi : 1 1 — > t and P^ is a one parameter semigroup due to homogenity of idi); similarly define 
Qra.lv an( j qW _ pij-g^iy^ note that by defintion of a and 6, 











■a 









t.t+l/n r i,t+l/n 



(y n W) (ii) 



for t being a multiple of T/n. Secondly it is intuitively clear, and an easy exercise to show, that 
P™/ = P t v _ s resp. Q"f = Q^ s for s, t points on the dissection D n . Since (a (rT 1 , ■) , 6 (n"\ •)) 
converges to the path id2 : t i — ► (t, t) in the sense of <j8j) (that is, £™ (t) = a and z n (t) = 

b tj, q = p = 1), we know that 

kv.w (0>2/o; (a («-~\ •) , & (n -1 , ■))) = y n ~> V = ^v,w (0, j/o;id 2 ) as n -> oo 
in ■Ioq.jq T j norm. Using the identity ([TT]) one recovers the "classic Lie-splitting" 

Qi/n ° Pf/raJ (l/o) -»• yt as n co for every t g [0, T] 

where jy is the ODE solution of (flO|) . Moreover, the convergence holds in .^.[qt] norm and by 
interpolation even in stronger (1 + e)-variation norm for every e > 0. 

Note that no rough path theory is needed and everything follows from continuouity in the sense 
of (|SJ) (with q = 1) which can be established by elementary computations; see |FV09[ Chapter 3]. 
Let us remark that using Young integration theory one can push this method to driving signals of 
finite p— variation for p ^ 2. Paths of Brownian regularity or worse are outside the scope of Young 
theory but one can use rough path results. 

2.2 Splitting RDEs 

Motivated by the above example we can ask for splitting for an RDE with drift of the form 

dy t = V (y t ) d£ + W (y t ) dz t , y (0) = y G R e , (12) 



5 



where z eC*"™* ([0, T] , (R<*)) and ^eC 1 "™ 1 ([0, T] , G^ (R d )) ; the path (£, z) has trivially "com- 
plementary Young regularity" and if V € Lip 7 (R e , K e ) and W = (Wi) C Lip'" 1 (R e , R e ) for 
7 > 1,7 > p then there exists a unique solution ttv.(w) (0, j/o; (£j z )) to (fT2|). For later use we 
show 

Lemma 1 Letz £C p ' var ([0, T] , G^ , £ 6 G«-" ar ([0, T] , R). If we define £ A (t) = £(a(A,i)) E 

C«-™ r ([0, T] , M), z A (t) = z(b (A, t)) e G p "™ r ([0, T] , G^ (R d )) then 



sup 1 1 z A 1 1 + sup 

A>0 ^ A>0 



< OO 



do (z A ,z) 



-> as A -> 0. 

oo;[0,T] 



Proof. First note that the variation norm is invariant under reparametrisation which implies the 
first statement. For the second statement it is sufficient to show pointwise converge (by intcr- 
palation pointwise convergence in combination with uniform variation bounds implies convergence 
in supremum norm). However, z and £ are, by assumption, both continuous paths which gives 
pointwise convergence. ■ 
Define 

:=t(a(n~\t)) and z? := z (b (n _1 , t)) . 

Similar to the ODE example, define the solution operator |p"' t ^,0 < s < t < Tj mapping points 

in R e to R e as P"' t V (x) :— Try (s,x;£ n ) t and the operators P v ,Q n,w and Q w . It remains to show 
that 

P£ v = and Q' s l f = QZ for s, t e D* 

(in contrast with the ODE example, P^ and are now two parameter semigroups due to the 
time-inhomogenity of £ and z). Since 

GM fad) 

is a geodesic space, there exists a sequence of paths 
(concatenations of geodesies on the sequence of dissections D m ), (z m ) m C G 1_var ([0, T] ,R rf ) with 
S[ p ] (z™t) = z s,t for s, t € L> m , such that 

sup||S [p] (z m )|| < oo, 

m 

do {S[p] ( z ™) 7 z ) — ► as m — > oo, 

with SVpi as in ([9]). Now define z n,m (i) := z m (a (n _1 , i)) and note that z m and z n,m have bounded 
1-variation (z m by construction, z n,m because the variation norm is invariant under reparametrisa- 
tion). Hence, for m,n fix we deal with an ODE as in the example above and therefore 

tt w (s, x; z"' m ) M = ir w (s, x; z m ) st 

for s,i 6 D n . Keeping n fixed and letting m — > oo, the LHS converges to 7i"vp (s, a;; z™) s t by lemma 
|l| and Lyons' limit theorem and the RHS to ttw ( s , x ', z ) s t ; we can conclude Q"' t = for 
s.t £ D n . A similar argument shows P"' ( V = P^ t for s,t S I?™. We finish the argument in 
the same way as in the previous example: solutions of ftvJW) ( s i x 'i {£" \ z ™)) converge uniformly to 
ttv,(w) ( s i x 'i (C; z ))- O n neighbouring points e D n , ttv.(w) (£ n , z ™)) s t can be identified as 



Qn,W -r> n -V 
St ° P s.t 



G 



Hence, for every t € [0,T] 



\t/n\-l 

n:Split 

Vt ■= 



II [Qfc/n,(fc+l)/n oP fc/n,(fc+l)/n (j/o) ^(ytQ (0, ^05 (f, z)) t »S Tl 00 a.S. 



fc=0 

Moreover, sup„ |?/™ ;Spllt | p _ var .r T ] < 00 which implies by interpolation convergence in (p + e)- 
variation norm of y™' S P ht for every e > 0. 

Remark 2 Similarly, one shows convergence of a splitting scheme, running the semigroups in the 
different order P^oQ^. Note also, that we restrict ourselves in this article to Lie splitting schemes 
but the methods can be easily modified to include Strang splitting (see \MQ 02l for the difference 
between Lie- and Strang splitting schemes) by using an appropiate modification of the time change. 
Further, we just deal with equidistant partitions. Numeruous variations of all this are possible ( as 
long as one can show convergence in a rough path topology of the approximating sequence) and such 
modifications are of great importance for rates of convergence - a topic which we hope to address in 
future work. 

3 Rough partial differential equations 

One could hope for similar results for SPDEs in Stratonovich form, 

d 

du = F (t, x, u, Du, D 2 u) dt + ^ A fc (t, x, u, Du) o dB k (13) 



fc=i 



u (0, x) = uq (x) on 



Indeed, splitting for equation (fT3")) has been treated in |KG03| for the case of F being a linear 
operator. To give meaning to (|13[) for nonlinear F one can introduce the concept of (rough) 
viscosity solutions (cf. jLS98al ILS98U ILS98bl ILSOOal ILSOObj and jCFO09j or [FUlO]). Let us 
informally discuss the idea before we give the precise definition in section 3) a real- valued, bounded 
and continuous function w on [0,T] x R e is called a solution if it is the uniform limit (locally on 
compacts) of (standard) viscosity solutions (u n ) of the equations 

d 

du" = F (t,x,u n ,Du n ,D 2 u n ) dC +J2A k (t,x,u,Du)dz n '\ 



i=l 



u n (0, x) = u a (x) on (0, T] x 



where (z n ) C C°° ([0,T] ,R d ) and (£ n ) C C°° ([0,T] ,M) are sequences of smooth driving signals, 
converging to a weak geometric rough path (£, z) (see ©)• Formally we write 

du = F (t, x, u, Du, D 2 uj d£ + A (t, x, u, Du) dz. 

(a solutions is then also a solution in the sense of jLS98bj when z is an enhanced Brownian motion). 
In the case when the rough path z is an enhanced Brownian motion and £ (t) = t, this gives a 
natural concept of Stratonovich solutions to (1131) . However, to apply the methods outlined in the 
sections above to derive a splitting method, care has to be taken: firstly, in the RDE case "any 
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sequence" (£",z™) of smooth paths converging to (£,z) gives rise to a solution, but in the RPDE 
case, if £ t < 0, one can not expect to treat even the simple second order equation 

du n = D 2 u n dC (14) 

since it is no longer degenerate elliptic (£ t < amounts to running the heat equation backwards in 
time). Secondly, assume we want to use Lie-splitting on dyadic partitions, i.e. approximate £ (t) = t 
with £ n (i) = a (n _1 ,ij_. This introduces a discontinuous time-dependence (£ does not exists on 
points of the partition ]Jj in equation (|14p . Such time-discontinuities are in general difficult to handle 
in a viscosity setting. Thirdly, one has to show continuous dependence of the solution^ of (fT5|) , not 
only on z but continuous dependence on (£, z) in a rough path sense. 

The first point is dealt with by characterizing the class of legit approximations £™ , leading to 
the path space Cg~ ml,+ ([0, T] , R), described in section |3~T1 Section 15721 deals with nonlinear PDEs 
with a discontinuity-in-time introduced by and section U gives the precise definitions of rough 
viscosity solutions and stability, contains the main theorem and examples of stable RPDEs. 

3.1 The space C 1_var ' + ([0, T] , M) 

As pointed out above, we have to avoid to fall outside the scope of (degenerate) elliptic PDEs. 
Using the notation C 0,1— var ([0, T] , R) for the closure of the space of smooth paths in variation 
norm (C^ l - ll -™ r ([0,T] ,R)) we recall that 

Wq 1 ([0,T],R) = jx : [0,T] ->R, By G L 1 ([0,T],R) s.t. x(t) = y (u) du 
= {x : [0,T] — >• R, x absolutely continuous, x (0) = 0} 

= {x : [0, T] — >■ R, x i C°°, x(0) = 0} Ml ~ var 
ee C°' 1 ~ uar ([0, T] , R) ^ Co 1 -™ 1 - ([0, T] , R) . 

Definition 3 C* 1,+ ([0,T] ,R) = G Cq 1 ([0,T] ,R) : ^ T =T,f > q|. 

Note that the paths a (A, •) and b (A, •) are not elements of Cg' + ([0, T] , R) , but elements of its 

closure, Cq :+ ([0,T] ,R) , in sup-norm. Working with Cq' + enables us in section below to give a 
short proof of existence, uniqueness and stability of a solution to PDEs of the type d t u = F£ t for 

paths £ G Cq' + '°° ([0,T],R). 

Proposition 4 Denote C^ var ' + ([0,T] ,R) = C 1,+ ' °° ([0,T] ,R). T/ien 
C i-« ([0!T])R) = |^ e C o ([0,T],R) :£ T -T and 3£ E L 1 ([0,T] ,R> ) , 

3a G CQ"" ar ([0, T] , R>o) , a increasing, a = a.s. and £ t = a* + / t; u du 

Jo 

7 One could avoid discontinuous time-dependence by restricting the class of splitting schemes (i.e. the class with 

C C 1 ). However, nearly all popular schemes (Strang, Lie, etc.) would then not be covered. 
8 The results in ICFO09| and [FO10 do not cover this due the time-discontinuity of the approximating sequence 

(<e). 
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and Cl~ var ' + ([0, T] , R) £ C^ var ([0, T] ,R). 

Proof. C: Let (£ e ) be a Cauchy sequence wrt. 1.^. Since (Co ([0,T] ,R) , j.j^) is complete, £ e 
converges uniformly to some £ G Co ([0, T] , R). This £ is monotone (not necessarily strict) increasing 
and hence |CI i— va r-[o t] < 00 ( recau that £ T = T). Every function of finite 1-variation is Lebesgue- 
a.e. differentiable and has a representation of the form 

£ t = a t + / £, u du 
J[o,t] 



where a is a function of 1-variation with a — Lebesgue-a.e. Now £ s . s+h > 0, for every h > 0; 
s G [0, 1). Hence we have a SyS+ h > — jg +h £ u d« and sending h — >• shows together with a = a.s. 
that a is monotone increasing and this implies £ u > Lebesgue-a.e. 

D: F (t) := £ t defines a continuous distribution function on [0, T] and let X be a random variable 
with distribution F. For e > denote by F e the distribution function of the random variable X + eN 
where N is a standard normal, independent of X. Clearly, X + sN — > X a.s. as e — >• and so 
the i^converge pointwise. By the lemma below this implies uniform convergence of F e to F. It 
remains to show that £f := (i) is C 1 but this follows from 

F E (t) = f F{t-u) dF eN (u) 
Jo 

where F E jy is the distribution function of eN. m 

Lemma 5 Let (f ri ) ri>0 C Co ([0,T] , R) ,f v (1) = l,each f v increasing (not necessarily strictly) and 
assume f r > -> / G C ([0,T] ,R) pointwise as n -> 0. T/ien, j/ 7 * - /l^.^ ->• as ->■ 0. 

Proof. Given e > we can choose a n G N big enough s.t. 

'G)-'(^)|<s 

for every « G {0, 1, . . . , n} . Now choose r? small enough such that 



n ) \n 



<2 



for all i G {0, 1, . . . , n) . This implies \f v (x) — f (x)\ < e since every x is an element of (at least 
one) interval , ^] and by monotonicity and using above estimates 

f- W </"(£)< + 1 + |s /(.) + .. 

Similarly 

/"(a:) >/(*)+£ 
and so (a;) - / (x)| < e for all x G [0, T]. ■ 

Remark 6 Concerning the choice of notation C L ~ t,ar ' + , note that the space of paths of finite 1- 
variation Cg _ " ar zs given as the closure of C 1 -paths with uniformly bounded 1-variation. Since 
paths in Cq' + have 1-variation bounded by T, the notation Cq var - + ([0,T] ,R) seems natural. 
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Remark 7 The paths a (A, •) and b (A, •) converge in (1 + e)-variation to the path id\ : t i — >• t for 
every e > and therefore also uniformly (but not in 1-variation!). 

Remark 8 G'l~ var ' + ([0,T] ,R) is not a linear space but a convex subset of CQ~ var ([0, T] , R). 

Remark 9 Despite the restriction of C^ var,+ ([0, T] , R) to paths with £ (T) = T which is con- 
venient in the proofs, one can handle PDEs with general increasing processes by rescaling; e.g. 
replace £ &j/ f (t) := £ (*) S C 1 ~ war ' + ([0, T] , R) and write dw = F (t, x, u, Du, D 2 u) d£ = 

F (t, x, u, Du, D 2 u) d\ with F := F^p- . 

3.2 PDEs with discontinuous time-dependence 

This section extends the notion of viscosity solutions to equations of the form 
du = F (t, x, u, Du, D 2 u) d£ (t) , u (0, x) = uq (x) , 

with F a continuous function and £ € CQ~ var,+ ([0, T] , R). In the appendix we show that this 
solution concept coincides with the notion of generalized viscosity solutions (going back to |Ish85j ) 
whenever the latter exists. In view of applications in sections section 0] and O and to keep tech- 
nicalities down, we focus on time-independent F. A proof for time-dependent F is given in the 
appendix. 

Proposition 10 Let (<f ) £ C ([0,T] ,R) converge uniformly to some £ G C^ var;+ ([0,T],M) 
as e — > 0. Assume {v £ ) e C BUG ([0,T] x R™,R) are locally uniformly bounded viscosity solutions 

of . £ 

d t v e = F £ (x, v s ,Dv e , £>V) £ t , v e (0, x) = v (x) . (15) 

with F £ : R e xRxR e x§ e — !• R, § e denoting the space of symmetric (e x e)-matrices, a continuous and 
degenerate elliptic function. Further, assume that F e converges locally uniformly to a continuous, 
degenerate elliptic function F and that a comparison result holds for dt — F E = and dt — F = 0. 
Then there exists a v such that 

v e v locally uniformly as e — > 0. 

Further, v does not depend on the choice of the sequence approximating £ and we also write v = 
to emphasize the dependence on £ and say that v solves 

dv — F (x, v, Dv, D 2 v) d£ t , v (0, x) = vo (x) . 

Prepare the proof with 

Lemma 11 Let £ G Cq' + ([0, T] , R), 

F : [0, T] x R™ x R x R n x S™ -> R 
and let F (t, x, r,p, X) = F (t),r,x,p,X) for (t,x,r,p,X) G [0, T] x R n x R x R™ x S™. Then 

1. if u G Bf/C ([0, T] x R n ) is a sub- (resp. super) solution of dt — F£ = 0, w(0, .) = uq (.) i/ien 
u> (f, x) := u , x) is a sm&- (resp. super) solution of dt — F = 0, uj (0, .) = uq (.). 



10 



2. if w G BUC ([0, T] x R") is a sub- (resp. super) solution of d t — F — 0, w (0, .) = wo (■) then 
u (t, x) := w (£ t , x) is a sub- (resp. super) solution of dt — F£, u (0, .) = wq (.). 

Proof. 1. Let ip G C 1 ' 2 ([0,T] x R") and assume w(t,x) — tp(t,x) attains a local maximum at 
(i,x) e [0,T] x R n . Then 

w (t, x) -cp (t, x)=u (C 1 (t) ,x) - ip (i, x)=u (C 1 (t) ,x)-ip (C 1 (i) , x) 

where ip (i, x) :— cp (£ t ~, x) . Using that u is a subsolution gives 

dt(p\q^ < f^'C 1 (t) ^(r 1 (t),x),D^^D^\^^ e c _ 1 

= F (x, r 1 (*) , w (*, *) , D<p\ u , D\\ l& ) 
= F (x, t, w (t, x) , Dtp^D 2 ^.) 

where we used that £></?|,-i . = £ and D 2 <£>L-i ~= -D 2 <£>|£ Since dt£>L-i = c?t<^|t a Clt- 1 an d 

S£i ' S£ ) ' St; 'St 

£ > it follows that 

dM lx <F(xXw{i,x),D^\ l& ,D 2 V \ i& ) . 

Now the same argument as above when u is a supersolution. 

2. Let G C 1 ' 2 ([0, T] x R") and assume u (t, x) — <p (t, x) attains a local maximum at (t, x) G 
[0,T] x R™. Then 

u (t, x) - 95 (i, x) = to ((5 (f) , x) - <p (i, x) = w (£ (i) , x) - (£ (f) , x) 

where </? (t, x) := y> (i) , x) . Using that w is a subsolution gives 

= F (x, i, u (i, x) , Dip\ it& , D 2 ip\ i £ ) 
where we used that Dip\ f -i -= D<p\$ % and D 2 ip\ f -i .= D 2 ^!^. Since £ > and d t (p\^ x = d t tp\t,x 

(r 1 )' ie t = &vk*(£ (*)) it; follows that 

d t <p\ lx < F(x,t,u(t,x) ,D ( p\ {j£ ,D 2 cp\ u )t(i) . 

Now the same argument as above when w is a supersolution. ■ 

Proof of Proposition HOi Set w e (t, x) := v e (t) , xY by lemma [TTT 

w £ is a solution of dt — F s £ = iff w e is a solution of dt — F £ = 0. 

Let 

w := limsup *w e and w :— liminf *w £ 

and note that 

F £ (x, r,p, X) — > i 7, (x, r,p, X) locally uniformly. 
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Standard viscosity theory tells us that w and w are sub- resp. supersolutions of dt — F = 0. Using 
the method of semi-relaxed limits (by definition, w > w and the reversed inequality follows from 
comparison which holds by assumption) conclude that w (t, x) := w (t, x) = w (t, x) . Further, using 
a Dini-type argument, for every compact set K C R n , 

\vf — ui|oo-[o T i xK = sup \w e (t, x) — w (t, x)\ — ► as e — >• 0. 

' te[0,T], iGK 

Now define 

v(t,x) := w(£(t),x) . 
and we get the claimed convergence v £ —> v. ■ 

4 Splitting RPDEs 

4.1 Rough viscosity solutions and stability 

Solutions of RDEs can be defined as limit points of ODEs. Similarly one can define solutions of 
rough partial differential equations as limit points of PDE solutions. 

Definition 12 Let z eC °^ mr ([0, T] , G® (R d )) , £ G C^ var;+ ([0, T] , R) and denote (£, z) G C °' p -™ r ([0, T] , (R d+r 
t/ie Young pairing (given canoncially via Young integration) . Further, let (£ £ ,z £ ) C Cq' + ([0,T] ,R)x 
Cq""" 1 " ([0, T] , R d ) , converge to (£, z) m i/ie sense o/ (HI) and assume the PDE 

d 

du £ = F (x, u £ ,Du £ , D 2 u £ ) df + ^ A fe (i, ac, u £ ,Du £ ) dz £;fe on (0, T] xl", 

fc=i 

u(0,ar) = u {x) e BUC(m e ,R) 

has a unique solution u £ € BUC ([0, T] x R e ;R) for every e > 0. We call every limit point u of 
(u £ ) (in BUC topology) a solution of the RPDE 

du = F (x, u, Du,D 2 u) d£ + A(i, x, u, Du) dz on (0, T] x R™, u (0, x) — u$ (x) . (16) 

If additionally, the limit is unique, does not depend on the choice of the approximating sequence 
(£ £ ,z £ ) and the map 

(£, z) G Cl~ var ' + ([0, T] ,R) x c°' p "™ r ([0, T] , G M (M d )) i — y u G BC/C ([0, T] x R e , R) 

is continuous then we say that the RPDE is stable in a rough path sense and we also write 
u = u^ ,z (or u — u z when t; (t) = t) to emphasize dependence on the rough path (£, z). 

4.2 The main theorem 

We are now able to formulate our main theorem. The proof is an easy consquence of the results in 
the previous sections. In section 0] we show that the assumptions are satisfied for a large class of 
RPDEs. 
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Theorem 13 Let z eC^ p ~ var ([0, T] , (R d )) ,£ £ C* +4 ~ t ' ar ([0, T] , R) and assume u e BC/C is 
ifte unique solution of the stable ( in the sense of definition \12\) RPDE, 

du = F (x, u, Du, D 2 u) d£ + A (t, x, u, Du) dz on (0, T] x R e , u (0, as) = u (x) e Bt/C. (17) 

Assume further that also the two (R)PDEs given by setting either F = or A = in < \17fy are 
stable. Denote {P s ,t, < u < T} and {Q S) t, < s < t < T} £/ie solution operators 

P M : (R e , R) — s> £>[/C ([0, T] x R e , R) , 931 — ► u, 

Q st : BJ7C*(R e ,R) -> BL/C([0,T] x M e ,R), 0i — > w, 



with 



and set 



Then 



dv = F (x, v, Dv, D 2 v) d£, v (0, x) — ip (x) , 
dw — A (t, x, w, Dw) dz, w (s, x) = <f> (x) . 

u n;Spte ^ ^ . = Y[ [Q in -x t{i+1)n -i o P in _ li(i+X)n _i] (uq (x)) . 
i=0 

u n\Spht _^ u i oca Hy uniformly as n — > 00. 



Proof. Define z" = z(&(n _ \i)), f(t) = ^(a( n -\t)). By lemma [TJ (£",z n ) -+n in the 

sense of (p} and by stability, the solutions u n of 

dw" = F (x, u n , Du n , D 2 u n ) dC + A (t, x, u n , Du 11 ) dz" on (0, T] x R", u (0, x) = u (x) 

converge to u, the solution of (|16p . Now for each given n one can identify on points of the dissec- 
tion {k—, i = 0, . . . ,n} the solutions of u ™; s P ht with u™ and by assumptions u n converges locally 
uniformly to u. ■ 

4.3 Examples of stable RPDEs 

This section shows stability in a rough path sense for a large class of RPDEs. Splitting results 
then follow readily by theorem 1131 Throughout this section z is a geometric p— rough path, i.e. 
z eC °" r ([0,T] ,Gp (R d )) , P > 1. 

Proposition 14 Lei 

L(x,r,p,X) = Tt\a(x) t x] +b(x) -p + f (x,r) 
K k {t,x,r,p) = (p ■ a k {t,x)) + rv k {t,x) + gk{t,x) 

with A (x) = a (x) a T (x) G S e and a : R e — > R exe , b (x) : R e — > R e bounded, Lipschitz continuous 
in x. Also assume that f : R e x R — > R is continuous, bounded whenever r remains bounded, and 
with a lower Lipschitz bound, i.e. 3C < s.t. 

f(x,r)-f(x,s) >C(r-s) for all r > s, x£R e . 
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and that the coefficients of A = (Ai, . . . , A^) , thats is o,v and <?, have Lip 1 -regularity for 7 > p + 2. 
Then the RPDE 

du = L (x, u, Du, D 2 u) dt + A (t, x, u, Du) dz, u (0, •) = Uo (18) 

is stable in a rough path sense and has a unique solution u z G BUC ([0, T] x R e ,R). 

Proof. We use the same technique of "rough semi-relaxed limits" as in [FO10]; for details of the 
transform of u £ to v £ we refer to |FO10| : the key remark being that u e is a solution of 

du £ = L (x, u £ , Du £ , D 2 u £ ) di £ + A (i, x, u £ ,Du £ ) dz £ , u (0, •) = u Q (•) (19) 

iff 

v £ (t, x) := {^y 1 (t, u £ (t, iP £ (t, x)),x)+a £ (t, x) ; 

is a solution of 

dv £ = L £ (x, i> £ ,Dv £ , D 2 v £ ) d£ e , v £ (0, •) = u (■) . 

Here L £ is a linear operator with coefficients determined by the characteristics of the PDE dw = 
A (t, x, w £ , Dw £ ) dz £ and tp £ , 4> 6 , a a are ODE flows converging to RDE flows <fi z ,a z which depend 
on z (but not the approximating sequence of (z £ )). Further a comparison principle applies to 
dt — L £ = 0. The assumptions of proposition [10] are then fulfilled, L £ L locally uniformly and 
using the method of semirelaxed limits, 

v £ — > i) locally uniformly. 

Unwrapping the transformation, that is, setting 

u z (t, x) := <f (t, v (t, (ry 1 (t, x)) - a z (t, (^-J" 1 (t, xj) ) , (20) 

finishes the proof since stability of the RPDE follows directly from the representation (f20|) . ■ 
Proposition 15 Let 

F(x,r,p,X)= inf \ Tr A(x\a) T X + b (x; a) ■ p + f (x, r; a) \ 

where A, a and b, f satisfy the assumption of proposition \14\ uniformly with respect to a G A and 
v = (f k)fc=i C Lip 1 (E n ,K n ) 7 > p + 2. Then the RPDE 

du = F (x, u, Du 7 D 2 uj dt + Du ■ v (x) dz, u (0, •) = uq 

has a unique solution u z G BUC {[0,T] x W 1 ) and is stable in a rough path sense. 

Proof. Similar to the proof above. ■ 
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5 Applications to stochastic PDEs 



The typical applications to SPDEs are path-by-path, i.e. by taking z to be a realization of a 
continuous semi- martingale Y and its stochastic area, say Y (us) = (Y,A); the most prominent 
example being Brownian motion and Levy's area. Taking the linear case as an example, the 
stability result of proposition [TJ] allows to identify 

du = L (t, x, u, Du, D 2 u) dt + A (t, x, u, Du) dz, u (0, •) = uq 

with z — Y (uS) as Stratonovich solution to the SPDE 

du = L (t, x, u, Du, D 2 u) dt + A (t, x, u, Du) o dY, u (0, •) = uo- 

Indeed, under the stated assumptions, the Wong-Zakai approximations, in which Y is replaced by 
its piecewise linear approximation, based on some mesh {0, — T\, the approximate solution 
will converge (locally uniformly on [0, T] x R™ and in probability, say) to the solution of 

du = L (t, x, u, Du, D 2 u) dt + A (t, x, u, Du) dY, u (0, •) = uq, 

as constructed in proposition 1141 In view of well-known Wong-Zakai approximation results for 
SPDEs, ranging from |BF95l[Twa95] to [GM041 IGS06] . the rough PDE solution is then identified as 
Stratonovich solution. (At least for L uniformly elliptic: the (Stratonovich) integral interpretations 
can break down in degenerate situations; as example, consider non-differcntiable initial data uo 
and the (one-dimensional) random transport equation du = u x o dB with explicit "Stratonovich" 
solution uq (x + B t ). A similar situation occurs for the classical transport equation ii = u x , of 
course.) Motivated by this, if u z is a RPDE solution of 

du = F (t, x, u, Du, D 2 uj dt + A (t, x, u, Du) dz, u (0, •) = uq, 
then we call u z with z = Y (w) as Stratonovich solution and write 

du = F (t, x, u, Du, D 2 u) dt + A (t, x, u, Du) o dY, u (0, •) = uo- 

The following example was suggested in [LS98c and carefully worked out in [BM02. BM07 . 

Example 16 (Pathwise stochastic control) Consider 

dX = b (X; a)dt + W {X; a)odB + V (X) o dB, 

where b,W,V are (collections of) sufficiently nice vector fields (with b,W dependent on a suitable 
control a = a (t) G A, applied at time t) and B,B are multi- dimensional (independent) Brownian 
motions. Defin$\ 



v (x, t; B) = inf I 



g{X*> t )+ / f(X^,a s )ds 



B 



9 Remark that any optimal control a (■) here will depend on knowledge of the entire path of B. Such anticipative 
control problems and their link to classical stochastic control problems were discussed early on by Davis and Burnstein 



DB92;. 



15 



where X x,t denotes the solution process to the above SDE started at X (t) = x. Then, at least by a 
formal computation, 

dv + inf [b (x, a) Dv + L a v + f [x, a)] dt + Dv ■ v (x) o dB = 

aeA 

with terminal data v (•, T) = g, and L a = ^ Wf in Hormander form. Setting u (x, t) = v (x, T — t) 
turns this into the initial value (Cauchy) problem, 

du = inf [b (x, a) Du + L a u + f (x, a)] dt + Du ■ V (x) o dBT-- 

aeA 

with initial data it (•,()) = g; and hence of a form which is covered by theorem \15\ f Moreover, the 
rough driving signal in proposition \15\ is taken as zt := By-t (w) where B (uj) is a fixed Brownian 
motion). 

Using theorem 1131 we immediately get a splitting result: 

Example 17 (Splitting HJB-equations) Let B be standard d— dimensional Brownian motion. 
Then the SPDE 

du = inf \ Tr a (x; a) a (x; a) T D 2 u + b (x; a) ■ Du + / (x; a) \ dt + (Du ■ v (x)) o dB, 

aeA I I i J 

u(u>;0,x) = u (x), (21) 

has a unique a unique solution u if o~(x,a) : R 6 x A —> R exe and b(x,a) : R e x A — > M. e are 
Lipschitz continuous in x, uniformly in a £ A, v = (v\,...,Vd) C Lip 7 (R e ; R e ) with 7 > 4. 
Denote {P u ,u > 0} the solution overato^\ of 

du = inf \ Tr a (x; a) a (x; a) T D 2 u + b (x;a) ■ Du + f (x;a)\ dt, (22) 
-yeA L L J J 

i.e. PiUo (.) = u (t, .) , and {Q s ,t, < s < t < T} the solution operator given by Q^.t^ (.) = ip (jr-v (s, .; B) t ) , 
with TT-y (s, x; B) t the SDE solution of 

dy=-V{y)odB,y s =xeR e . (23) 

L*/"J-1 



(t,x):= J] [Q 

i/n,i/n+l/n Pl/n] («0 {x)) U (t, x) as n -S> OO 



i=0 



and the convergence also holds locally uniformly. 

Thus, equation (|2ip can be approximated by solutions of a standard HJB equation (f2"TT) and by 
solutions of the RDE (1231) (for numerical schemes for HJB, see FS06a]~). 



10 Note that in the case £ = t one can use a one-parameter semigroup since = 2 on the time interval on which 
the approximation evolve. 

11 Apriori the leftmost two terms would have to be Q 1 t/ n j ,t Pt- It/nj n- However, the claimed convergence follows 
from Lyons' limit theorem. 



1G 



Example 18 (Linear SPDEs, Filtering) Let L and A be as in proposition^^ Then there exists 
a unique solution to 

d 

du = L (x, u, Du, D 2 u) dt + J2^k (t, x, u, Du) odB k ,u (0, .) = u (.) 

k=l 

Denote by {P u ,0 < u < T} the solution operator 

Lp i — s- v with v solution of dv = L (x, v, Dv, D 2 v) dt, v (0, •) = tp (•) 
and by {Q s .t,0 < s < t} the solution operator 

ip i — > y with y solution of dy — A (t, x, u, Du) o dB, v (0, •) = tp (•) 

Then, 

Lt/nJ-l 

u n;Spht ^ . = J| [Qi/nj/n+x/n o P 1/n ] (u (x)) ->u(t,x) as 12 -X OO 
i=0 

and the convergence also holds locally uniformly. 

A Appendix: Time-dependent F 

To deal with time-dependent F we need the additional assumption of uniform bounds on the 
derivatives of the approximating sequence £ e . 



< oo, converge uniformly to some 

oo;[0,T] 



Proposition 19 Let (<f) e C Cq' + ([0,T] ,R) , sup £ f 

£ E Cl~ var;+ ([0,T] ,M) as e -x 0. ylsswxie (V) £ C BUC ([0,T] xR e ,I) are Zoca% uniformly 
bounded viscosity solutions of 

d t v e = F £ (t, x, v £ ,Dv £ , D 2 v £ ) v £ (0, x) = v (x) . 

with F £ : [0,T]xl e xIxI e x § e — x K a continuous and degenerate elliptic function. Further, 
assume that F £ converges locally uniformly to a continuous, degenerate elliptic function F and that 
a comparison result holds for dt — F £ = and dt — F = 0. Then there exists a v such that 

v £ —X v locally uniformly as e —X 0. 

Further, v does not depend on the choice of the sequence approximating £ and we also write v = 
to emphasize the dependence on £ and say that v solves 

dv = F (t, x, v, Dv, D 2 v) d£ t , v (0, x) = vo (x) . 

Proof. Set w E (t,x) := v E ^(£ e ) 1 (t) ,x^j and divide [0, T] into intervals on which £ is strictly 

increasing resp. constant, i.e. = sx < ti < S2 < • • • < t n = T, £ strictly increasing on [sj,ij], 
constant on [£(,S{+i]. By lemma [TT1 on intervals [sj,ii] 

v £ is a solution of 9< — F£ — iff w £ is a solution of <9f — = 
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where F £ (t,x,r,p,X) = F ((£ £ ) 1 (t) , x, r,p, x). Let 

W := limsup *w e and u; := liminf *w e 

e 6 

and note that on intervals [sj,£j] 

F £ (i, x, r,p,X) -4 F (i) ,x,r,p,X) =: F (t,x,r,p,X) locally uniformly. 

Standard viscosity theory tells us that on intervals [sj,ij] , W and w are sub- resp. supersolutions 
of dt — F = 0. Using the method of semi- relaxed limits (by definition, w > w and the reversed 
inequality follows from comparison) conclude that w := w = w and that for every compact set 
K C R™ (by using a Dini-type argument), 

\w e — w]^, t A xK = sup \w 6 (t,x) - w (t,x)\ -> as e ->• 0. 



Now define 

w (ew,x), r 1 (*) < < < r 1 (*o 

w{Z(ti),x), C 1 (U) < t < C 1 (s i+ i) 



v (t, x) :- 



We get the claimed convergence v e v on intervalls [£ (si)>£ 1 (i,)l . However, on [£ (ii)>£ 1 ( s i 
w e is by definition viscosity solution of dt — F^ = with initial condition v 6 (ij) , .) and 



F£ — > £ locally uniformly since sup t 



£ (t) 



is uniformly bounded in s by assumption. Hence 
the standard stability result of viscosity theory applies and v e converges locally uniformly on 
[C _ ( s i+i)l against the constant-in-time function v e (£«) > ■) > the only solution to 

9t = with initial condition v E {U) >■)■ This proves the claimed convergence. Further, note 
that v is given as the unique viscosity solution of dt — F = 0, hence every other sequence approxi- 
mating £ will lead to the same limit. ■ 

The proof of the main theorem (theorem I13[) and applications to examples adapt now in a 
straightforward way to time-dependent F. 



B Appendix: Generalized viscosity solutions 

Section 13.21 and appendix extend the notion of viscosity solutions to equations of the form 

du = F (t, x, u, Du, D 2 u) d£ (i) , u (0, x) = u Q (x) . (24) 

with £ G C 1 -" ar;+ ([0,T] ,R) . Generalizations of viscosity solutions go back to |Ish85j . |LP87j and 
for the parabolic case |Nun92j . Let us recall the definition given in [Nu n92j . 

Condition 20 F (-,x, r,p, X) G L 1 ((0,T),R) for all (x,r,p,X) 6i e xlxl e x§ e and F is 
continuous on R e x R x R e x S e for almost all t £ (0, T) 

Condition 21 F (t, x, -,p, X) is nondecreasing on R for all t G (0, T) and for all (x,p,X) G 
R e xR e x S e . 
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Definition 22 Let F satisfy conditions \20\ and \21l A locally bounded uniformly upper semicontin- 
uous function u G BUC ([0,T] x R e ) is called a generalized subsolution of 

du = F (t, x, u, Du, D 2 u) , u (0, x) = uq (x) (25) 

if for any (i,x) S [0,T] xW,be L 1 ([0,T] ,R),0 G C 2 (M e , K) , G :[0,r]xM e xMxl e xS e -)M 
continuous and degenerate parabolic, such that 

u (t, x) + b (r) dr — <p (t, x) attains a local maximum at (t, x) , 
Jo 

and 

b(t) + G(t,x,r,p,X) < F (t,x,r,p,X) for a.e. t G B s (i) and 

for all (x, r,p, X) G Bg (x, u (t, x) , D(f>\^ F) 2 ^ ^ for some 8 > 

it follows that 

b (£) + G (f , u (t, i) , £>V| t - 4 ) < 0. 

A locally bounded uniformly lower semicontinuous function is called a supersolution if the above 
estimates hold when one replaces maximum by minimum and reverses the inequality sign. 

Note that equation (|25[) is covered by this definition. However, it is quite cumbersome to derive 
existence, comparison and stability results in this very general setting and in the case of interest to 
us, the time-discontinuouity only appears multiplicatively. 

Proposition 23 Under the assumptions of vrovosition [W\ and additionaly £ 6 W 1 ' 1 the function 
u = is a viscosity solution of 

du = F (t, x, u, Du, D 2 uj d£ t , u (0, x) — uq (x) 

in the sense of definition \22{ 

Proof. We partition [0, T] into < si < ti < ■ ■ ■ < s n < t n < T such that £ is increasing on 
[sj, ti\, constant on [ti, Sj+i]. Say u (t, x) + J b (r) dr — <j> (t, x) attains a local maximum at (t, x). 
If t G [si,tj] by construction u(t,x) = w(t,£ t ) with w a viscosity subsolution of dt — F = 0, 
F (t,r,x,p,X) — F (£ _1 (t) ,r,x,p, X) , hence also a generalized subsolution and using that £ is 
invertible on [sj,tj] one sees by a change of variable that also u is a generalized subsolution on 
[ Si ,ti] of d t -F = 0. 

If t G [ti,Si + i], then £ is constant, hence £ = a.s. and so F (t, x,r,p, X) £ (t) — for a.e. 
t G Bs (t) and u is a generalized subsolution on that interval. This shows that u is a generalized 
subsolution and the same argument shows that u is a generalized supersolution. ■ 
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